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Abstract 

The problem of explaining the acceleration of the expansion of the universe 
and the observational and theoretical difficulties associated with dark matter 
and dark energy are discussed. The possibility that Einstein gravity does 
not correctly describe the large-scale structure of the universe is considered 
and an alternative gravity theory is proposed as a possible resolution to the 
problems. 
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1 Introduction 

The surprising observational discovery that the expansion of the universe is ac- 
celerating PU El IS] has led to an increasing theoretical effort to understand this 
phenomenon. Although interpreting the data by postulating a non-zero, positive 
cosmological constant is the simplest way to understand the data, it is not satis- 
factory, because it leads to the two serious problems of why the estimates from the 
standard model and quantum field theory predict preposterously large values of the 
cosmological constant 4J, and the coincidence of dark energy dominance today. 

If we simply postulate a repulsive force in the universe associated with a charge 
density, then we might expect that this force could be responsible for generating 
the acceleration of the universe. However, for a homogeneous and isotropic universe 
the net charge density would be zero, although for a finite range force with a small 
mass there will exist a non-zero charge density The effect of a Maxwell-type 
force would be to lower or raise the total energy, leaving the form of the Friedmann 
equation unchanged. Thus, we would still have to invoke exotic forms of energy with 
an equation of state, p = wp, where w is negative and violates the positive energy 
theorems. For a non-zero cosmological constant w — — 1. 
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In addition to the dark energy problem, we are still confronted with the puzzle of 
dark matter. Any observational detection of a dark matter candidate has eluded us 
and the fits to galaxy halos using dark matter models are based on several parameters 
depending on the size of the galaxy being fitted. The possibility that an alternative 
gravity theory can produce a satisfactory description of galaxy dynamics with little 
or no dark matter should be studied El El CD! • 

Challenging experimental results are often the precursors of a shifting of scientific 
paradigms. We must now entertain the prospect that the discovery of the mechanism 
driving the acceleration of the universe can profoundly change our description of the 
universe. Recent observational data for supernovae (SNe la) [2] have produced new 
conclusive evidence that the universe went from a decelerating phase at z ~ 0.5 to an 
accelerating phase and that the data are consistent with the concordance model with 
Q m ~ 0.3, Da ~ 0.7, where Q m and fl\ denote the ratios of dark matter density and 
dark energy density to the critical density p c = 3H 2 /8tcG. Moreover, the data give 
for the equation of state parameter for dark energy, wd = —1.02 ± 0.15 consistent 
with a cosmological constant value wd = — 1 and dwo/dz = 0. 

Given the uneasy tension existing between observational evidence for the accel- 
eration of the universe and the mystery of what constitutes dark matter and dark 
energy, we shall consider the question of whether Einstein's general relativity (GR) 
is correct for the large scale structure of the universe. It agrees well with local 
solar system experimental tests and for the data obtained for observations of the 
binary pulsar PSR 1913+16 [TT]. However, this does not preclude the possibility 
of a breakdown of the conventional Einstein equations for the large-scale structure 
of the universe. The standard GR cosmological model agrees well with the abun- 
dances of light elements from big bang nucleosynthesis (BBN), and the evolution 
of the spectrum of primordial density fluctuations, yielding the observed spectrum 
of temperature anisotropies in the cosmic microwave background (CMB). Also, the 
age of the universe and the power spectrum of large-scale structure agree reasonably 
well with the standard cosmological model. However, it could be that additional 
repulsive gravitational effects from an alternative gravity theory could agree with all 
of the results in the early universe and yet lead to significant effects in the present 
universe accounting for its acceleration. 

A fundamental change in the predictions of the observational data will pre- 
sumably only come about from a non-trivial alteration of the mathematical and 
geometrical formalism that constitutes GR. Recent modifications of Einstein's grav- 
itational theory and developments of cosmological models ^21 EH [13 HH| have 
led to alterations of the Friedmann equation at large cosmological scales. 

In the following, we shall consider the physically non-trivial extension of GR 
called the nonsymmetric gravitational theory (NGT) ^7j. This theory was exten- 
sively studied over a period of years, and a version of the theory was discovered 
that is free of linearized weak field inconsistencies such as ghost poles, tachyons, 
and other instabilities jTHj. The skew sector corresponds to a massive Kalb- 
Ramond field, and is not to be identified with the electromagnetic field as was done 
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by Einstein [T2j in his unified field theory. The skew sector is treated as part of the 
total gravitational field and together with the symmetric part Qi^u) forms the total 
nonsymmetric g^ U) which generates a non-Riemannian geometry. 

A significant modification of the infrared gravitational field and the equations 
of cosmology can only be obtained by having additional degrees of freedom beyond 
the two degrees of freedom of Einstein's gravity theory. Such additional degrees of 
freedom are supplied by NGT, which yields six additional degrees of freedom in the 
skew sector g\p U ]- 

As we shall see in the following, and in further work j20], it is possible for 
NGT to describe the current data on the accelerating universe and the dark matter 
halos of galaxies, gravitational lensing and cluster behavior, as well as the standard 
observational results, without invoking the need for dominant, exotic dark matter 
and a dark energy identified with vacuum energy and the cosmo logical constant. 
These results are obtained while retaining the good agreement of Einstein's gravity 
theory with solar system tests, terrestrial gravitational experiments and the binary 
pulsar PSR 1913+16. 

2 NGT Action and Field Equations 

The nonsymmetric g^ u and are defined by [T7J [TH] : 

fiV = 5W) + 9[M (1) 



and 
where 



— + rj^„], (2) 
g( M u) = -{g^v + gun), #[h = -(g^u - g V p). (3) 

The contravariant tensor g^ u is defined in terms of the equation 

g^gau = g u »gu« = 5",. (4) 

The NGT action is given by 

^ngt — S + Sm, (5) 

where 

S = J d 4 x[g^R; u (W) - 2Av^ - ^ 2 g^ M ], (6) 

and Sm is the matter action satisfying the relation 

1 ( 5S M \ _ _1 



r =g\5g> lv J 2 
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Here, we have chosen units c = 1, = y/—gg^ u , g = Det(<7^), A is the cosmological 
constant, /x is a mass associated with the skew field gW|. Moreover, is the 
nonsymmetric energy-momentum tensor and R* (W) is the tensor 



r; v (w) = r pv (W)-^w,w u , (8) 

where R pu (W) is the NGT contracted curvature tensor 

RAW) = W% 3 - \{W^ U + W^) - W*,Wfr + W&WZ, (9) 
defined in terms of the unconstrained nonsymmetric connection: 

= rj, - (io) 

where 

H/„ = 

2 



^ = ^(^A-^)- (11) 



Eq. lfTUJ) leads to the result 

r„ = rf M] = o. (12) 

The contracted tensor -R^W 7 ) can be written as 

R^(W)=R^(T) + ^W M , (13) 
where ^ 

-R^(r) = r^,/3 - 2 (^(8),^ + ^(vp),ij) ~ r^r^ + r^r^. (14) 

A variation of the action S^gt yields the field equations in the presence of matter 
sources 

G;„{W) + Ag,„ + = SitGT^, (15) 
g [H , = -\g {lMa) W a , (16) 

+l%srwfa + \{z m W p 5» - t>WWp%) = 0. (17) 
Here, we have G*^(W) = R*^(W) - \g^K*(W), where 11* (W) = g» v R*^(W), and 

Su* = -^ 2 (#[h + \g„ v g [ap] g\p*] + g [ap] g^g P u). (18) 

The generalized Bianchi identities 

[g av G pv (T) + g«G vp (Y)\ a + g»\ p G» v = 0, (19) 
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give rise to the matter response equations 

g^\ v + g m T v ^ v + ( w + g pVttx - g mp )T^ = 0. (20) 

It has been proved that the present version of NGT described above does not 
possess any ghost poles or tachyons in the linear weak field approximation ^Hj, either 
as an expansion about Minkowski spacetime or about a generic GR background. 
This cures the inconsistencies discovered by Damour, Deser and McCarthy in an 
earlier version of NGT |21j . In NGT there are three distinct possible metric tensors 
with three different local light cone structures. The definition of a spacelike surface 
is consequently dependent on the chosen coupling of matter to geometry and it is not 
possible to unambiguously apply a (3 + 1) decomposition of field variables in order 
to perform a Hamiltonian constraint analysis using the standard methods. Further 
studies of the non-linear and non-perturbative solutions of massive NGT and its 
Cauchy development have to be undertaken. In the following, we shall identify g^ u ) 
with the metric tensor of spacetime. 



3 Cosmological Solutions 

For the case of a spherically symmetric field, the form of g^ u in NGT is given by 
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where a, /3, 7 and w are functions of r and t. The tensor g^ v has the components: 



\ T^— 







P 2 +P 

/csce 

'/3 2 +/ 2 







/csce 

/3 2 +/ 2 
/3CSC 2 6> 

/3 2 +/ 2 





w*— 07 





7 



(22) 



We have 



^ = s in6[(a 1 -w 2 )((3 2 + f 2 )] 1 / 2 . 

For a comoving coordinate system, we obtain for the velocity vector 
which satisfies the normalization condition 

n J_ 



g^ v )U^u v 



u 



u' 



u 



= 1: 
0. 



(23) 
dx^/ds, 

(24) 



We set w = so that only the g[23] component of g\nv\ is different from zero. This 
corresponds to setting the magnetic monopole charge and static magnetic field in 
Maxwell's theory to zero, if we define 



*\pu] 



(25) 
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where e fMUap is the Levi-Civita tensor density and we identify g*^ and g*M k ] as the 
static "magnetic" and "electric" potentials, respectively, associated with the massive 
Kalb-Ramond potential field. 

The vector can be determined from 

= -^g^, (26) 

where sr^g^ 01 ^ = 5 u fl . For the skew symmetric field with w — 0, it follows from 
(|IB|I and that W 7 ^ = 0. 

The energy-momentum tensor for a fluid is 

T» v = (p + p)u"u v - pg^ + K.M , (27) 

where is a skew symmetric source tensor identified with the intrinsic spin or 

fluid vorticity. The fluid vorticity is defined by 

= U M + Q V l H> ( 28 ) 

where a M is the fluid's four-acceleration. We only consider couplings to U[^ v ] to avoid 
derivative couplings. A rotational action is given by 

S R = J d 4 xC R = J d 4 x K pg^e^u [a , p] , (29) 

where k is a coupling constant. By varying this action with respect to g^ u we get 

fiua/3 

KM = Kp —=u M . (30) 

We define 

Tfiv = 9,/39auT al3 , (31) 

and from (jlj) and (|27j). we get 

T = p — 3p + g [a/3] K^ =p-3p + 2fK, (32) 

where K^ = Kj sin 9. 

We can write the field equations (|15|) for = in the form 

1 1 

R^(T) - g^A + - -g^S = 8irG(T^ - -g^T), (33) 

where S = g^S^. The metric takes the canonical Gaussian form for comoving 
polar coordinates 

ds 2 = dt 2 - a{r,t)dr 2 - (3{r,t)[d9 2 + sin 2 9 <d<j) 2 }. (34) 



6 



We shall assume that ft{r,t) ^> f(r,t) and that a solution can be found by a sepa- 
ration of variables 



a(r, t) = h(r)R 2 (t), (3{r, t) = r 2 S 2 (t). 
From the field equations (|55jl (see Appendix A), we get 



R S _ 1 
R ~ S ~ 2 ' 



(35) 



(36) 



where R = dR/dt and Z is given by 



z= ftT_ 5ftft'f 2 _ aPT + Wff ff 3/'/ 



ft 3 



2ft 4 



2aft 3 



ft" 



ft 2 2ft 2 



, dff + 2/3'// 



2a/? 2 ' /?3 • (37) 

We shall assume that Z ~ which from (j3T?j) gives J?(t) ~ S'(i). This leads to a 
metric of the form 

ds 2 = dt 2 - R 2 {t) \h{r)dr 2 + r 2 {d8 2 + sin 



(3* 



We cannot impose exact homogeneity and isotropy on the skew sector , for this 
would lead to f(r,t) = w(r,t) = 0. 

We shall further simplify our calculations by assuming that the mass parameter 
/i ~ 0, that the cosmological constant A = and that we can neglect any effects 
due to the antisymmetric source tensor associated with vorticity of the matter 
fluid. 

With the assumption that ft ^> /, the equations of motion become (see, Ap- 
pendix A): 



2b(r) + R(t)R(t)+2R 2 (t)-R 2 (t)W(r,t) = 4irGR 2 (t)\p(r, t) - p(r, t)}, 



R(t)R(t) + l -R 2 (t)Y{r,t) = ^R 2 ( t )[p(r,t) + 3p(r,t)}, 



where 



2b(r) 



h'{r) 
rh 2 (r) 



(39) 
(40) 

(41) 



The functions W and Y are given by 



ur( oc'ft'f 2 ft"f 2 aftf 2 5ft' 2 f 2 aff 
W{r, t) = - — + -^—^ + 



2a 2 ft 3 aft 3 2aft 3 2aft A 
a'ff ff" Aff ft' , 3f 2 



2aft 2 



2a 2 ft 2 aft 2 aft 3 2a ft 2 ' 



(42) 



7 



Within our approximation scheme, W and Y can be expressed in the form 

W ( tu h ' f2 2/2 2k2f2 10/2 kf i h ' ff 

vv[r,t) -+- -+- 



/ l 2 r 5 jR 6 h r 6 R 6 r 4 R 6 fa&jp r 4_R5 2/i 2 r 4 i? 6 

//" sff + J/^ (44) 



M 4 i? 6 /ir 5 i? 6 2/ir 4 i? 6 ' 

1 ' ; r 4 i? 6 r 4 i? 6 2r 4 i? 4 r 4 i? 5 r 4 i? 4 ' 1 ; 

Eliminating R by adding (|3U|) and pUj). we get 

#(t) +&(r) = ^p(r,t)i2 2 (t) + Q(r, (46) 

where 

From (|4U|) we obtain 

= -^^)[p(r,t) + 3p(r,t)) + ~R{t)Y(r,t). (48) 

Let us consider an expansion of f(r,t) about a background fo(r,t): 

f(r,t) = f Q (r,t) + 5f(r,t) + .... (49) 

We shall identify the fluctuations Sf(r, t) about the background fo(r, t) as any matter 
content additional to the visible baryon matter of the universe, p m (r, t) = p(Sf(r, £)). 
Thus, p m replaces the exotic cold dark matter (CDM) of the standard cosmological 
model. We also have the expansions 

W = W + 6W + Y = Y + 6Y + ....; 

Q = Q + 5Q + .... (50) 
We define the matter density to be 

PM=Pb + Pm, (51) 

where p& denotes the baryon density. The background field / will describe the 
source of dark energy density and we shall consider slowly varying solutions of fo, 
which rise towards the present epoch with red shift z ~ and a solution of 5f that 
yields a p m = p(Sf) that decreases with increasing time as 1/R 3 . 
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We can now write in place of (}4*6j) and (}48|) : 

R 2 {t) +b(r) = ^p M (r,t)R 2 (t) + Q (r,t)R 2 (t), (52) 

and 

R(t) = -^R(t)\p M (r,t)+Zp M (r,t)] + ±R(t)Y (r,t), (53) 

where Q Q {r,t) = Q(f (r,t)) and Y (r,t) = Y(f (r,t)). 
We can write Eq. ljB^jl as 

H 2 + ^ = MI 2 , (54) 

where H = R/R, 

Q = Q M + Q fo , (55) 

and 

8nGp M Q 

M = ~!ip~' /0 = W (56) 

From Eg. (|41jl for h = 1, we obtain 6 = and f2 = 1 and 

H 2 = ^ P m + Qo, (57) 

which describes a spatially flat universe. The line element takes the approximate 
form of a flat, homogeneous and isotropic FRW universe 

ds 2 = dt 2 - R 2 (t)[dr 2 + r 2 (d6 2 + sin 2 9d<j> 2 )]. (58) 



4 Accelerating Expansion of the Universe 

Let us assume that there exists a solution of the field equations such that asymp- 
totically for large r and t: 

Q' o (r,t)~0, y o '(r,f)~0, (59) 

and 

p M (r, t) -> p M (t), p M (r, t) -> pm(*)- (60) 

Here, we choose the big bang to begin on a hypersurface S(r, t) with r = t = 0, so 
that asymptotically for large r and t we approach the present universe. 

It follows from that R > when F > ^G{.Pm + 3^a/)- If we assume that 
there exist solutions for Q (r,t) and lo(r, t), such that they are sufficiently small 
in the early universe, then we will retain the good agreement of Einstein gravity 
with the BBN era with p rad oc l/R 4 . As the universe expands beyond the BBN 
era at the temperatures, T ~ 60 kev-1 MeV, and p M oc l/R 3 , then we must seek 
solutions such that Q and Y begin to increase and reach slowly varying values with 



9 



f2° o ~ 0.7 and Q° M ~ 0.3, where Q° M and f2° denote the present values of Qm and 
Qf , respectively. Provided we can find solutions of the field equations that satisfy 
these conditions, then it should be possible to fit the combined supernovae, cluster 
and CMB data 0. 

We observe from ({H]) and that Qo and Y are functions of the behavior of 
R and / and their derivatives. If fo grows sufficiently with Rast increases, then Qo 
and Yq can dominate the matter contribution pu as the universe evolves towards 
the current epoch with pu « 0. A detailed solution of the field equations is required 
to determine the dynamical behavior of R, f, Q and Y. 

In the present epoch, p M m and gives 

^ = -±[4nGp M -Y ]. (61) 

We can now define an effective equation of state parameter for the universe: 

w eS = ^(A7rGp M -Y ). (62) 

We see that Y > AmGpu corresponds to the usual condition for acceleration W£> < 
— 1/3 with the dark energy equation of state pd = wdPd- 

We can explain the evolution of Hubble expansion acceleration within NGT, 
without violating the positive energy conditions for matter and radiation. We satisfy 
the strong energy condition for matter pu + %Pm > throughout the evolution of 
the universe and there is no need for a cosmological constant. The Qo an d Y 
contributions to the expansion of the universe increase at a slow rate up to values 
today with Y > AixGpM (pm ~ 0). The cosmological constant A = during this 
evolution. 

We must also guarantee that the influence of the Y and Q contributions do 
not conflict with galaxy formation, i.e. the additional NGT contributions to the 
Friedmann equation must not couple too strongly with the attractive gravitational 
effects predicted by the field equations in the galaxy formation epoch. This issue 
and the other required evolutionary effects of Q and Y must be determined by a 
numerical computation of the NGT field equations. 

5 Conclusions 

We have proposed that NGT may explain, as gravitational phenomena, the acceler- 
ating universe without unknown dark energy, and the observed flat rotation curves 
of galaxies without the undetected exotic dark matter [20] • We do expect that there 
is some dark matter in the universe in the form of dark baryons and neutrinos with 
non-vanishing mass. Such a theory can be falsified with data, whereas it is difficult 
to falsify the dark energy and dark matter hypotheses. 

NGT would be required to explain the formation of galaxy structure without 
CDM. We know that the standard ACDM model of structure formation and the 
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description of the CMB power spectrum is remarkably successful, so that solutions 
of the NGT equations must produce an equally successful description of the data. 
To investigate in detail whether NGT can succeed in producing a successful account 
of the large-scale structure of the universe will require solving the field equations. 
These are issues that require further investigation. 



6 Appendix A: The T-Connections and Field Equa- 
tions 

The NGT compatibility equation is given by 

g\u )V - g P u^ p x v - 9x P ^u = \^ pp) (9 P ug\ v - g w gx P - gxug[ PV ])W^, (63) 

where is determined from (}26|) . When w(r,t) = g[oi] (r,t) = 0, it follows that 
W/j, = and the compatibility equation reads 

g\u,r, - g P uT p Xv - g\ P V p nv = 0. (64) 

The non-vanishing components of the T-connections are: 
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rfi 3] = -Ifi2] sin 2 q= x -b sine, 

rfao] = - r [20] sm 2 = ~£>sin0, 
sin 9 / 1 
~2tV2' 



where A, £? and Z? are given by 

A = ln(/3 2 + / 2 ), 
P 2 + f 2 ' 

and 

n _ g/ - /g 
/? 2 + / 2 ' 

The NGT field equations in the presence of sources are given by 

R n (T) = — A - -\(A ) + 4B 2 } + + — - 

/VV dfa\ a/7 a 1 a . 1 2 af 2 
"UtV + ^l^J + 2^l^-^ + 2 A J +Aa -4 /i ^T^ 
= 47vGa(p-p + 2fK), 

R2 2 (T) = R, 3 (T)cosec 2 9 = 1 + f^"^' )' + (^f^) + V«) 
B(fA' + 2(3B) 8 /2fD-(3A\ 2fD-j3A 

D 1 Bf 2 

--(/A + 2/?£>) + A/3 + ~/^j 2 = 4nG[3(p -p + 2fK), 

Roo(T) = - l -A -\A 2 + AD 2 ) + ±A + *- f I - *- 

oov ; 2 8 V ; 4 7 2a V2 7 2a 



dt\2aj \2aJ 2a\2a 2 7 2 
R m (T)=8nGK [10] =0, 

12 



R [23] (T) = sin 6 
B 



J 8a \a 7 



4a 



_(2/ B -/M<)-£(M + 2^)(i + £ 



9t V 47 / 47 



-47rG/sin6'(p -p). 
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